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Slow dynamics of the contact process on complex networks 
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Abstract. The Contact Process has been studied on complex networks exhibiting dif- 
ferent kinds of quenched disorder. Numerical evidence is found for Griffiths phases and 
other rare region effects, in Erdos Renyi networks, leading rather generically to anoma- 
lously slow (algebraic, logarithmic,...) relaxation. More surprisingly, it turns out that 
Griffiths phases can also emerge in the absence of quenched disorder, as a consequence 
of sole topological heterogeneity in networks with finite topological dimension. In case 
of scale-free networks, exhibiting infinite topological dimension, slow dynamics can be 
observed on tree-like structures and a superimposed weight pattern. In the infinite size 
limit the correlated subspaces of vertices seem to cause a smeared phase transition. These 
results have a broad spectrum of implications for propagation phenomena and other dy- 
namical process on networks and are relevant for the analysis of both models and empir- 
ical data. 
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^ 1 Introduction 

£f) ■Nonequilibrium systems have been a central research topic of statistical mechanics [£T]— [3|] . As in equi- 
t— I librium they exhibit various phase transitions, critical phenomena and universality classes. Since 
^ 'they cannot be treated by thermodynamical formalism universality provides a guidance's in the zoo 
of models exhibiting diverging correlation length and scaling behavior [4]. Power-law scaling be- 
havior has been reported in many real and model systems for spontaneous occurrence the theory of 
self-organizing-criticality was proposed [5]. This, however requires the competition of a slow accu- 
mulation and fast dissipation mechanism. Diverging correlation length on the other hand is generated 
naturally in nonequilibrium systems by the currents or fields acting on them. Furthermore, quenched 
disorder may also cause extended regions in the parameter space with non-universal power-law dy- 
namics [44]. A fundamental dynamical system model is the Contact Process (CP) 00], in which sites 
can be either occupied (infected) or empty (susceptible). By changing the infection rate of the neigh- 
bors A/k, where k is the degree of the vertex, a continuous phase transition occurs at the A c critical 
point from inactive to active steady state. The inactive state, characterized by the order parameter, the 
density of infection (p) is zero. This is also called absorbing state, because no spontaneous activation 
of sites is allowed. 

In the new century the interest is shifting from models, defined on Euclidean, regular lattices to 
processes defined on general networks 10 |9[. Since the introduction of a simple model describing 
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the emergence of scaling in random networks fl 1 Ofl the study of complex networks is flourishing. 
Multidisciplinary applications involve, for example, the WWW, various biological, sociological and 
technological networks. These are, mostly scale-free (SF) networks, exhibiting P(k) ~ k~ 7 degree 



distribution of the nodes (for recent reviews see [8),[ 1 1 J). Other families of complex network models 
are those composed of a af-dimensional regular lattice and additional long edges 01211 . These arise in 
the context of conductive properties of linear polymers with cross-links that connect remote monomers 



1 13], in public traffic systems [14], in the case of nanowires n social [16] and phone networks 



1711 or to describe forest fires 111 811 among other examples. In these models a pair of nodes separated 
by the distance I of the base lattice are re-connected by a long edge with the asymptotic probability 
for large /: 

Pd)=pr s - (i) 



In the special case 5 = 0, edges exist with a length-independent probability, as in small world net- 
works, therefore these models are called Generalized Small World networks (GSW). The s = 2 case is 
important from network optimization point of view lfl4ll for example. If s > 2, they are characterized 
by a finite topological dimension d, i.e. N(l) ~ l d , where N(l) is the number of nodes within the 
graph distance / from a given node. Note that if the links vary quickly with respect the time scale 
of the dynamical process defined on these networks they realize Levy flights i20ll2lll . Usually these 
annealed networks exhibit different behavior than the quenched ones. 

Although many network models exhibit infinite topological dimension (d), simple mean-field ap- 



proximations cannot capture several important features 112 3142811 . Very recently it has been conjectured 
II29M3 ill that generic slow (power-law, or logarithmic) dynamics is observable only in networks with 
finite d. This claim is relevant in the light of recent developments of dynamical processes on complex 
networks such as the simple model of "working memory" 1I32I . brain dynamics ll33ll . social networks 
with heterogeneous communities 13411 . or to understand the slow relaxation in glassy systems [35]. 
Slow dynamics has been shown to originate from the bursty behavior of the agents connected by 
small world networks resulting in memory effects l42ll . On the other hand it can also be related to 
arbitrarily large (I < N), correlated rare-regions (RR), which possess long lifetime in the inactive 
jhase, above the pure critical point A® < A < A c . This can be understood by non-perturbative methods 



phase, ; 



More recently the possibility of power-law dynamics of CP has been investigated on different BA 



networks with y — 3 I28L I45H . Extensive simulations showed p(t) oc l/(fln(f)) density decay and 
p(A, t — > 00) oc \A — 1 1 steady state behavior with logarithmic corrections in agreement with the HMF 
approximations. On loop-less BA trees the epidemic propagation slows down and a nontrivial critical 
density decay emerges p(f, A c ) oc r° 5 . Additionally, when k dependent weighting was applied, which 
suppress hubs or make the network disassortatative GP-like dynamics was observed in the simula- 
tions. However, systematic finite scaling study revealed that these power-laws saturate in the N — > 00 
thermodynamic limit, suggesting smeared phase transitions known from Euclidean, regular systems 
if the correlated subspaces can undergo phase transitions themselves, when they are effectively above 
the lower critical dimension of the problem: dm? > d~ B44I1 . In this case, the dynamics of the locally 
ordered RR-s completely freezes, and they develop a truly static order parameter. Clearly in infinite 
dimensional networks such RR-s can be embedded as a percolation analysis confirmed this [45]. 



2 Optimal fluctuation theory 



The basic idea of the optimal fluctuation theory is that the long-time decay of the order parameter 
(p(t)) is dominated by the regions of size I, which are rare in general: P(l) oc exp(-c/), but can exhibit 
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exponentially long lifetimes t(1) oc exp(bl). In particular for the the density of infected sites of CP, 
RR-s provide the leading order contribution in the inactive phase: 



p(0~ J lP(t)exp(-t/T)dl, (2) 

which in the saddle point approximation results in p(t) ~ t~ c ' decay, with a non-universal, A de- 
pendent exponent 12911 . This is the so called Griffiths Phase (GP) H43L 14411 . bordered by a stretched 
exponential decay law at the critical point A® of the clean system, and by the dirty critical point 
A c > A®, where the evolution becomes logarithmically slow 

P (t) ~ Ht/toT & , (3) 

thanks to lim^^ b — 0. The a and other activated scaling exponents are described by the strong 
disorder universality fixed point in the renormalization sense. Extensive numerical simulations in 



1 3111 provided estimates for them in case of GSW-s. 



3 Griffiths phases in Generalized Small World network models 

GSW-s in one dimension have the intriguing feature that in the marginal case (s — 2) intrinsic prop- 
erties exhibit power-law behavior and the corresponding exponents vary continuously with the pref- 
actor B. First it was shown that the topological dimension of such networks (see Fig. [Tall depends on 
P JH-H. It has been claimed in a recent letter Q and in HEO] that if d(j3) is finite, GP-s and 
similar rare-region effects can also appear. Density decay simulations, started from fully active state 



(p(0) = 1), have been presented in 03111 . The effective decay exponents, defined as the local slope of 
p{t): 

\n[p(t)/p(Q] 
]n(t/f) 

saturate to A dependent constant values in the long time limit (Fig.fTbb. although logarithmic correc- 
tions arise. 
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Interestingly such GP-s have also been found in 3 -regular random networks, constructed as shown 
J49I1 . For larger values of s the long -ranged links are irrelevant and d(J3) = 1, while for s < 2 the 
topological dimension diverges and mean-field behavior emerges. 

4 Griffiths phases in Erdos Renyi network models 

The CP on Erdos Renyi (ER) graphs [51] with a quenched disordered infection rates (QCP) has been 
studied in [29-31]. A fraction q of the nodes (type-II) propagate infection with a reduced value Ar, 
with < r < 1, while the remaining fraction 1 - q (type-I nodes) take their "clean" value A. Pair 
mean-field approximations lead to the following critical threshold: 

^) = 7# L TT j -' (5) 

(k)-l l-q 

Type-I nodes experience a percolation transition, where the type I-to-type I average degree is 1, i.e. 
at q pe rc — 1 - (k)~ l - For q > q perc activity cannot be sustained: type-I clusters are finite and type-II 
ones do not propagate activity. Optimal fluctuation theory and simulations show that in this case the 
absorbing phase of QCP splits and a GP with power-law dynamics appears for A > A®. 

In fact in fragmented ER networks, where d is zero the intrinsic disorder of As is not necessary, 
but the topological heterogeneity itself is enough for the occurrence of GP-s as shown in Fig. [2b] 
Density decay simulations from fully active state of the CP on ER graphs with (k) — 0.15 up to sizes 
N = 10 6 confirm this. For high values of A one can see plateaus on the p(t) curves, similarly as 
reported in [27]. These are the consequence of the metastable local active domains. 

5 Slow dynamics in Barabasi-Albert network models 

After showing the rare-region effects and the possibility of GP-s on GSW and ER networks let us turn 
towards the results on the most common type of scale-free (SF) networks. CP have been simulated 
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and analyzed by mean-field methods on BA networks, in particular for loop-less and weighted cases 



as described in H45H . BA construction is a simple and fast way to generate SF networks, in contrast 
with other standard models, e.g [52]. The BA growth starts with a fully connected graph of small 
(iVo = 10) nodes. Following that at each generation step s, a new vertex with m edges is added to the 
network and connected to an existing vertex s' of degree k s > with probability YI S ^ S > = k s >/ 2. S "<. S kf. 
This process is iterated until reaching the desired network size N. The resulting network has a SF 
degree distribution P(k) k~ 3 and for m — 1 we obtain a BA tree (BAT) topology, while for the 
looped case m = 3 is used. Binary (non-weighted) BA networks can be transformed into weighted 
ones by assigning to every edge connecting vertexes i and j a symmetric weight cotj. In 14511 two 
different network topology dependent weight assignment strategy was introduced in order to slow 
down and localize epidemics. 

(i) Weighted BA tree I (WBAT-I): Multiplicative weights, suppressing the infection capability of 
highly connected nodes 

ojjj = a> (kikjy v , (6) 

where a>o is an arbitrary scale and v is a characteristic exponent with v > 0. This can model internal 
limitations of hubs, like the sub-linear Heap's law ll54Tl . 

(ii) Weighted BA tree II (WBAT-II): Disassortatative weighting scheme according to the age of 
nodes in the network construction 

H ~ J\ x n , 
"ij = —fi-, (7) 

where the node numbers i and j correspond to the time step when they were connected the network. 
Since the degree of nodes decreases as kj oc (N/i) 1 ^ 2 during this process, this selection with x > 
favors connection between unlike nodes and suppresses interactions between similar ones. 

The presence of these weights affects the dynamics of the CP. Thus, the rate at which a healthy 
vertex i becomes ill on contact with an infected (active) vertex j is proportional to Ao>y, therefore the 
epidemic can in principle become trapped in isolated connected subsets. Density decay simulations 



1 45], started from fully active state show GP like regions both in the WBAT-I and WBAT-II cases. For 
WBAT-II the scaling appears for t > 10 4 MCs in the region A > 9.5, even in case of networks with 
x = 3 as shown on Fig. [3a] The effective decay exponents, using t/f = 8, saturate to A dependent 
constant values in the long time limit (see inset in Fig. Oat. However, by increasing N ata given A the 
decay curves saturate asymptotically, suggesting a smeared phase transition as shown on Fig. [3b] 



6 Discussion and Conclusions 

I have overviewed the effects of quenched heterogeneity on the dynamics of the Contact Process in 
different network models. In finite topological dimensional cases slow dynamics (power-law, loga- 
rithmic, stretched exponential ... etc) and Griffiths Phases can be observed, which is very important 
for understanding spreading-like phenomena of real world networks. In these cases topological dis- 
order of the network can result is GPs. Similarly, in models where the roles of space and time are 
exchanged "temporal Griffiths phases" have been found in various systems near their phase transition 
points H El. 

On the other hand on graphs of infinite topological dimensions like in weighted SF trees, the 
network heterogeneity can again cause slow dynamics of the CP. However, in the thermodynamic 
limit the power-laws saturate and the locally active, high dimensional subspaces cause smeared phase 



transitions [45]. Very recently slow dynamics has also been reported in case of the Susceptible 



Infected-Susceptible model defined on the (SF) flower model II27I1 . weighted BA trees H28I1 and on ER 



graphs with exponential weights H60I1 . Spectral analysis of the quenched mean-field theory [26, 28] 
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Figure 3: |(a)| Density decay as a function of time in CP defined on WBAT-II trees, with exponent 
x = 3. Network size: N = 4 x 10 5 . Different curves correspond to: A = 9.4, 9.5, 9.6, 9.65, 9.7, 9.8, 9.9 
(from bottom to top). Inset: the corresponding local slopes |(b)l Density decay as a function of time in 
CP on a weighted WBAT-II trees with exponent x — 2 at A = 6.8. Different curves correspond to sizes 
N = 4 x 10 3 , 4 x 10 4 , 2 x 10 5 , 10 6 , 10 7 (from bottom to top). Inset: the corresponding local slopes. 



and simulations confirm epidemic localization effects in such models. Further studies of different 
models and networks are currently under way. 
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